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EFFECT  OF  GRADIENTS  IN  MULTI-AXIAL  STRESS  STATES  ON 


RESIDUAL  STRESS  MEASUREMENTS  WITH  X-RAYS 


by 

I.  C.  Noyan* 


ABSTRACT 

The  assumption,  that  stress  components  in  the  direction  of  the  surface 
normal  are  negligible  in  traditional  residual  stress  determination  methods 
by  x-rays,  has  been  recently  disproved.  In  this  paper  we  investigate  the 
effect  of  normal  stresses  on  the  accuracy  of  these  traditional  methods. 

It  is  shown  that  appreciable  error  can  exist  in  surface  stresses  determined 
by  such  methods,  if  normal  stresses  are  present.  New  procedures  are  proposed 
to  minimize  these  errors. 


^Graduate  student  and  Research  Assistant,  Dept,  of  Materials  Science  and 
Engineering,  Ihe  Technological  Institute,  Northwestern  University,  Evanston 
Illinois  60201 


INTRODUCTION 


The  determination  of  surface  residual  stresses  in  poly -crystal line 

(1  2) 

materials  by  X-ray  methods  is  well  established.  *  '  Traditional  methods 

assume  a  bi-axial  stress  state  that  is  uniform  in  the  surface  layers 

penetrated  by  the  x-ray  beam.  This  assumption  is  based  on  the  conception 

that  this  penetration  depth  is  too  low  to  be  affected  by  the  stresses  in 

the  third  dimension.  These  methods  predict  a  linear  variation  of  inter- 

planar  spacing  HdH  with  "sin3^" ,  with  a  slope  that  is  proportional  to  the 

stress  in  the  measurement  direction  (S^);  here  f  is  the  angle  between  the 

normal  to  the  diffracting  planes  (L^)  and  the  surface  normal  of  the  speci- 

men  (S^)»  as  shown  in  figure  1.  If  the  components  of  the  (assumed)  bi-axial 

stress  tensor  exhibit  steep  gradients  in  the  volume  sampled  by  the  x-ray 

beam,  curvature  occurs  in  the  "d"  vs  sin3^  plot.  The  slope  of  a  least-square 

fit  to  such  data  yields  the  stress  in  the  (S*)  direction,  averaged  over  the 

(3)  9 

depth  of  penetration. 

The  assumption  in  these  traditional  methods,  that  stress  components 

in  the  direction  of  the  surface  normal  (S~)  are  negligible  in  the  volume 

J  (4  5) 

sampled  by  the  x-ray  beam  has  been  disproved  in  several  recent  studies.  *  ' 

Stress  components  of  appreciable  magnitude  in  this  direction  have  been 

detected  in  the  surface  layers  of  materials,  created  In  various  ways. 

Recent  theory,  as  well  as  experiments,  show  that  if  stresses  (^3^23^ 

the  direction  of  the  surface  normal  (S^)  are  present,  "splitting”  of  the 

"d"  vs  sin2^  data  results;  that  is  "d"  vs  sin2T|r  plots  have  opposite  curvature 

for  negative  and  positive  if*  Analysis  of  such  data  is  described  in  detail 

by  Dc5lle,^  (see  also  reference  4  for  an  example  of  the  use  of  this  analysis). 

If,  however,  only  normal  stress  (^3)  is  present  in  the  direction  of  the 

surface  normal  and  shear  components  are  negligible,  curvature  occurs  in  the 

"d"  vs  sin2\|f  data  and  there  is  no  "^-splitting".  The  degree  of  such  curvature 

depends  on  the  steepness  of  the  gradient  in  (Since  i*  a  stress 

normal  to  the  surface,  it  is  zero,  by  definition  at  the  surface.  Thus  it 

has  to  exhibit  a  gradient  in  the  surface  layers.)  Dolle,  James  and  Cohen ^ 

have  shown  how  to  estimate  the  gradient  using  this  curvature.  It  will  be 

shown  in  this  paper  that,  even  for  very  small  curvature  in  the  "d"  vs  sin2* 

data,  analysis  with  the  bi-axial  assumption  causes  an  appreciable  error  in 

the  calculated  surface  stress.  Methods  are  proposed  that  provide  for  the 

determination  of  the  complete  stress  tensor  in  the  surface  layers  in  the 


1 


absence  of  {-splitting,  i.e.  when  <jjj>  ^3  are  ne8ligihle-  These  methods 
are  tested  here  by  computer  simulation* 


THEORY 

The  orthogonal  coordinate  systems  used  in  the  following  discussion 
are  shown  in  figure  1*  The  specimen  axes  are  defined  such  that  S ^  and 
are  in  the  surface  of  the  sample*  The  laboratory  system,  in  which  diffrac¬ 
tion  is  sampled,  is  defined  such  that  is  in  the  direction  of  the  normal 
to  the  family  of  planes  (hkl)  whose  "d"  spacings  are  being  measured  by 
x-rays*  L2  is  in  the  plane  defined  by  £^and  and  makes  an  angle  $  with 

S2-  In  what  follows,  primed  tensor  quantities  refer  to  the  laboratory 
system  L.  and  unprimed  quantities  to  the  sample  system  ^ ,  following  the 
convention  established  by  Dolle.  ' 

If  the  unstressed  lattice  spacing  "d0"  is  known,  strains  in  the 
direction  may  be  obtained  from: 


d-d 

(s')  »  -fy  -  ’ 

w33>$,+  d0 


(1) 


axes : 


This  strain  may  be  written  in  terms  of  strains  in  the  S.  system  of 

(6) 


(c^)^  ^  ”  e1^cos2$sin2ijr  +  e12sin2$sinsi|r  +  e,9sina$sinai|r  + 


’22 


e33cosa^r  +  e13cos§sin2^r  +  e23sin$sin2*  .  (2) 

Expressing  the  strains  on  the  right  side  of  equation  (2)  in  terms 
of  stresses  we  obtain: 


<‘33>M 


■  ^  {®nC0s2*  +  tf12sin2$  +  a22 8ina«  -  o33J  sin2*  + 

Or)  ff33~  I<*11+  ®22+  »33>  +  TT  {*13C0S*  +  ^ 23s in2+ 


Here  E  is  Young's  modulus  and  v  is  Poisson’s  ration.  If  a  stress  tensor  of 


the  form: 


'ij 


o  o\ 
o  022  O 


(4) 


-  2 


o 


o 


exists  in  the  volume  sampled  by  the  x-ray  beam,  equation  (3)  will  yield  the 
classical  residual  stress  equation: 

(*33 >#.*  “  ^  {'  11C08S*  +  *22s*“a*}sin3V  -  |(On  +  o22)  .  (5) 

Since  Cj^cos3$  +  *  Oj(the  stress  in  the  direction  S^),  equation  (5) 

becomes : 

^  sin3*  •  E  ({7ll  +  ff22)*  <6> 

It  is  seen  from  equations  (6)  and  (1)  that,  the  slope  of  MdM  vs  sin3i|r 

is  proportional  only  to  0,  when  a  bi-axial  stress  state  is  assumed.  If  a 

9 

tri-axial  stress  tensor  of  the  form: 


is  present  in  the  coordinate  system^,  equation  (3)  becomes: 

(<?33>M  “  ^  kncos2*  +  «r22sin3#  -  ^^sin3*  + 


E  ^11  +  ff22  +  ^°33^  ’ 


<<?33>«,^  “  ^  *  <a33>^  Slna*  +(1E^<ff33\ 


i  ^°ll  +  °22  +  ^33^-  * 


where  carats  "<  >”  imply  averages  over  the  depth  of  penetration  of  the 
X-ray  beam. 

From  equations  (8b)  and  (1)  it  can  be  seen  that,  when  a  stress  tensor 

such  as  equation  (7)  exists  in  the  surface  layers,  the  slope  of^a  linear 

least-squares  fit  to  "d"  vs  sin3^  will  be  proportional  where 

is  the  normal  stress  averaged  over  the  penetration  depth  (eqn.  26). 

A 

Thus  when  (0^3)  18  present  in  the  penetrated  volume,  there  will  be  an  error 
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r 


T 


in  the  stress  determined  from  the  slope  by  the  conventional  (bi-axial) 
analysis,  equal  in  magnitude  to  (0^3) • 

It  is  possible  to  obtain  both  (a and  using  equations  (3b)  and  (1). 
There  are  a  total  of  three  unknowns:  Stresses  and  $» 

which  is  the  angle  between  the  principal  axis  and  the  measurement  direction 
S|«  A  simple  method  of  determining  the  stresses  (0^3)  and  is  presented 
below. 

If  data  for  Md  vs  sin%  is  obtained  at  two  $  angles;  and  $A+  90°, 
where  is  arbitrary,  the  slopes  of  "d11  vs  sinaf,  +^qo  will  be: 


m*  ”  {^n0033* 


A  +  °228i“a*A  •  <ff33>] 


“  Of)  {®»  ”  <*33>}  •  (9b) 

A 

V  »°  ■  (?r)  t°uco,a  <V  +  <V  900  >  -  <»?3>}  <9c> 

Summing  Eqns,  (9a)  and  (?c)  we  obtain: 

W^°“^^11+522 " 2^33^  ’  <10a) 


0r;  (oll+  a22>  “  (ifc)  {\+  ^+90°}  +  2<®33>  •  (10b) 

The  intercept  I  of  ($^3)  in  Eqn,  8a  vs  sin3^  at  i)*GP  does  not  depend  on  f, 

1  "  (Te“)  <'a32>  "  I  ^ff33^  “  E  (C11  +  °22)  *  (11) 

The  constants  and  *  £  are  "X-ray  elastic  constants'*  and  may  be 

obtained  from  the  literature  for  a  given  reflection ,  hkl ,  and  material,  or 

-  4  - 


measured  following  the  procedures  outlined  in  reference  7.  (Since  only 
slopes  of  load  vs  "d"  are  used  in  these  procedures,  effect  of  gradients 
are  eliminated.)  Equations  (11)  and  (10b)  are  then  solved  simultaneously 
for  <0^3)  and  *22)*  With  ^33)*  the  stresses  in  and  +^<f 

directions,  a*  and  a  0,  are  obtained  from  equations  (9b)  and  ^9d) 
respectively. 


The  analysis  described  here  yields  the  stress  normal  to  the  surface; 

/A 

(oj j)»  and  two  stresses  in  the  surface  plane;  and  <*$A+ggo.  It  is  also 
possible  to  obtain  the  first  stress  invariant  Ij*  .  If  the  principal 

stresses  in  the  surface  5 22  an<*  C^e  an8^e  •  these  form  with  the  measure¬ 

ment  direction  are  required,  ”d"  vs  stnsf  values  must  be  obtained  for  a 
third  §  angle,  *A+  a,  where  a  is  arbitrary  ( a  ±  90°)  The  slope  of  the  "d" 

vs  sin2*  line,  M.  ,^will  be; 

$A-Wf 


1+V) 

\-H*  "  T  tff  llCOs3  <V  a)  +  a22Sin8  <V  “>  ‘  (a33>I  *  <12> 

.A  . 

Since  (a3y  is  known,  3^,  a 22  and  $  may  be  obtained  by  simultaneously 
solving  equations  (9a),  (9c)  and  (12). 

Special  Case 

Some  important  surface  treatments  such  as  shot  peenlng  produce  residual 
stresses  which  are  isotropic  in  the  surface.  The  stress  tensor  is  then  of 
the  form: 


(13),  where  6<<  Oy  • 


For  such  cases,  the  analysis  is  simplified  and  the  time  of  measurement 
is  reduced.  Substituting  equation  (13)  into  equation  (3): 

**33^  "(tt)Ki  “  ^33^}  Sln  *  +C~E>^33> 


“  E  {^33^*  +  29 ll} 


Since  the  slope  and  intercept  of  "(€33)^  vs  sinrT  may  be  solved  for 
°11  and  ^°33^»  one  arbltrary  *  tilt  is  sufficient  to  obtain  the  complete 
stress  tensor  in  the  surface  layers  sampled  by  the  x-ray  beam.  In  other 
words,  the  usual  measurement  of  f,d"  vs  sin2^  suffices,  but  the  slope  and 
the  intercept  at  f  ■  0 P  must  be  employed. 


Effect  of  Curvature 

The  term  "023"*  used  in  tbe  Previous  equations,  is  a  stress  normal 
to  the  surface.  It  is,  of  course,  exactly  zero  at  the  surface.  Therefore 
°33  Can  exist  only  as  a  Sradient  in  the  near  surface  layers.  It  is  well 
known  that,  when  such  stress  gradients  are  present,  curvature  occurs  in 
the  ndn  vs  sinsi(f  plots  and,  as  mentioned  above,  a  method  has  been 

proposed  to  use  this  curvature  to  obtain  the  gradient.  It  has  also 


proposed  to  use  this  curvature  to  obtain  the  gradient.  It  has  also 

been  suggested  that  the  slope  of  a  linear  least-squares  fit  will  yield  an 

(3) 

"average"  value  of  the  stress  state  over  the  depth  of  penetration. v  Si; 


Since 


the  x-ray  beam  penetrates  to  different  depths  at  each  ^-tilt,  a  different 
portion  of  the  existing  gradient  is  sampled  with  each  ^-tilt,  causing  the 
curvature  in  "d"  vs  sin2^.  Analysis  of  such  data  by  linear  least-squares 
may  cause  appreciable  error  in  the  slope  and  intercept  of  the  line  fitted 
to  the  data.  This  error  is  now  examined. 


law*: 


We  assume  that  the  behavior  of  stress,  0 with  depth  follows  a  power 
o  (z)  “  0^(2  •  0)  +  a^z^j  >  (15) 


Where  *  0)  is  the  stress  value  at  the  surface,  a_,  n^  are  constants 

(over  the  depth  of  penetration)  and  2  is  the  distance  coordinate  along  the 
normal  to  the  surface  (measured  into  the  material). 

The  average  stress  obtained  by  x-rays  penetrating  to  a  depth  "z"  is : ^ 


*For  example,  in  case  of  autof rettage,  (shrink- fitting  of  cylinders 
on  solid  shafts),  n  -  2  for  tangential  and  radial  stresses. 


6 


(17a) 


where  r  is  given  by 


.(7) 


sin* 6  -  sin£t 
2^sinGcosy 


for  ^  tilts  around  the  Q  axis  and 


(17b) 


sinQcosjfr 

Tt  2n 

for  tilts  around  an  axis  parallel  to  the  diffractometer  plane.  Here  ^  is 
the  linear  absorbtion  coefficient,  0  the  Bragg  angle  of  diffraction  and  if 
is  defined  in  figure  1.  Variation  of  and  with  are  shown  in  figure  2. 


Substituting  equation  (15)  into  equation  (16)  and  integrating: 


<*i 


CTij(2“0)  +  KijT  u 


(18) 


Where  is  a  constant,  a  is  the  exponent  of  the  stress  variation  (Eqn.  15) 
and  t  is  given  by  equations  17. 


Assume,  for  simplicity,  that  a  stress  tensor  of  the  form  given  in 
equation  (13)  exists  in  the  near  surface  layers  and  that  the  stress  in  the 
surface  plane  is  not  dependent  on  depth.  For  c33>  equation  (18)  becomes: 

<033)^  «  K33Tn33  (19),  where  n33±0  . 

Substituting  equation  (19)  into  equation  (14)  and  re-arranging: 


d_ 

do 


S3>,  -(t) 


2  K,,Tn33 

Pj^sin2*-^  «u>‘* - 1 - {(1+v)  (l-sin2i|;)  -vT  (20) 


(8) 

The  equation  of  a  least-squares  line  is; 


Y  *  B  +  8.X  +  e  , 
O  1 


(21) 


where  Y:  the  dependent  variable 

X:  the  independent  variable 

e:  error 

g  g  :  regression  parameters, 
o  1 

Comparing  equations  (6),  (20)  and  (21),  we  observe  that,  if  a  least- 
squares  line  is  assumed,  (i.e.  assuming  only  a  bi-axial  stress  state),  the 


error  term  will  be: 


K  ^33 

e^  *  - - - ( (1+v)  (l-sinat)  -  v!  . 


(22) 


It  is  seen  from  equation  (22)  that: 
i)  e  ^  is  a  biased  error.  (It  is  systematic.) 

-  7  - 


ii)  For  a  fixed  depth  (fixed  t),  the  error  "e  •'  in  (c  )  increases  with 

♦  33  * 

increasing  exponent  of  the  stress  (n^)* 
iii)  For  a  fixed  n^  *  there  will  be  greater  bias  in  (€^3)^  measured  at 

low  if  (high  t)  tilts  than  for  (*33)^  measured  at  high  f  (low  r)  angles 

Analysis  of  data  corresponding  to  equation  (20)  by  conventional  methods 
ignores  the  error  components,  e^ ,  associated  with  each  point  and  forces  a 
linear  least-squares  solution  to  the  data.  Since  the  error  is  not  random, 
but  highly  biased,  such  a  procedure  may  yield  large  errors  in  the  stress 

calculated  from  the  slope,  which  depend  on  a^*  ^3  in  equation  (15)* 

• 

When  a  tri-axial  stress  state,  (Eqn.  14),  is  employed,  there  will 
still  be  some  error  due  to  this  curvature;  expanding  equation  (20): 

<«33>t*(ii^{ffll"K33Tn33}Sina,,r  E^"33  +  *11}  (23) 

A  least-squares  fit  to  the  data  (with  equation  14)  yields  the  equation: 

^*33^“  (“T^®ll‘^°33^sinS*  +(.  E  )^33^“  E  ^033^+2<7ll^  * 


Thus,  the  deviation  of  each  point  from  linearity  is: 

JJ  f 


(K^rn33-<c^>) 


■{ (l+v)  (1-sin2 i)r)  +  v}  , 


where  (033)  is  the  least-squares  average  of  (^33)^  values* 


★The  equation  for  is: 

A 


<a33>-an- 


^<*1  -X)(Yi-  Y) 


iSfV3^ 


where  X  *  sin2f 


Y  *  N  *  Number  of  f -tilts  and 


/N  \  / 
X  -f  Z  X.\/N 

l1”1  if 


We  see  that  the  normal  stress  obtained  from  the  slope,  (-33)*  is  a  double 
average;  first  averages  are  taken  by  x-rays  over  different  regions  of  the 
stress  gradient,  than  these  are  averaged  again  by  the  least-squares  fitting. 


Comparing  equations  (25)  and  (22)  we  observe  that,  while  the  absolute 
magnitude  of  bias  associated  with  a  given  (el-),  is  less  when  the  tri-axial 

Jj  y 

analysis  is  used,  some  error  will  remain,  causing  an  error  in  stresses  cal¬ 
culated  from  the  slope  and  intercept  of  the  "d"  vs  sin2^  data.  These  errors 
will  be  termed  ’’curvature  errors"  here.  It  is  seen  from  equation  (25)  that, 
for  high  f  angles  (low  penetration  depths),  the  bias  associated  with  (*33)^ 
will  be  smaller  than  for  low  ^  points.  Thus  a  least-squares  analysis  of 
the  high  t  range  will  produce  a  more  meaningful  fit  of  the  linear  model  to 
the  data,  in  the  presence  of  such  curvature. 


General  Case 

Up  to  this  point  it  has  been  assumed  that,  no  shear  stresses  in  the 
direction  of  the  surface  normal  exist*  This  assumption  is  not  necessary 
as  far  as  the  curvature  error  is  concerned.  When  a  general  stress  state 
is  present  in  the  near-surface  layers,  the  strains  in  the  (j^)^  ^  direction 
are  given  by  equation  3.  The  analysis  in  this  case  has  been  described  by 


Do  lie: 


(6) 


The  quantities  a1 ,  a2  are  formed; 


(W 


[cjjCos2*  +  c12  sin2$  +  u^sin2* 


* 


sin2f  +  ^  £  ^(033)^  "  g  (ffll  +  c22  +  (27a) 

as  *  "  (^){<013>tcos4  +  <c23\Sin$-  Sinl2*l  ’  (27b) 


Where  e^+  is  the  e33  measured  t  *  +t  and  $,  and  e^-  is  the  e^3  measured 
at  i  *  -^1  and  $ . 


It  is  seen  that,  aA  and  a2  will  have  curvature  when  plotted  against 
sin2i(r  and  sin|  2^ |  respectively,  since  both  contain  components  that  are  a 
function  of  depth  (a^j  ^) .  In  fact  ax (Eqn.  27a)  is  a  general  form  of 

Eqn.  (24).  Using  high  i  -points  to  form  the  arrays  a!  vs  sin2*  and  a2  vs 
sin  |2i^  |  will  minimize  the  curvature  error  associated  with  this  analysis. 
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If  gradients  are  present  in  all  components  of  a  general  stress  tensor, 
equations  (20)  -  (25)  will  have  additional  terms..  However  the  general  ideas 
discussed  above  will  still  be  valid  for  most  cases. 


DATA  SIMULATION  AND  ANALYSIS 

2 

Computer  simulation  of  "dn  vs  sin  if  was  used  to  obtain  information 

about  the  magnitudes  of  errors,  resulting  from  the  presence  of  (0^3)  and 

the  curvature  error.  A  computer  program  was  written  which  calculated  the 

strains  in  the  JL  coordinate  system  for  a  given  stress  tensor,  converted 

these  strains  into  strains  in  the  (L0)x  .  direction  and  formed  the  "d"  vs 

^3  if 

sin2^  array  for  a  given  $.  The  generated,  C33  *  vs  sin2^  were  then  analyzed 
by  the  methods  described  above. 


Three  different  stress- tensor  groups  were  investigated.  The  assumptions 
common  to  all  three  groups  are: 

1)  The  hypothetical  sample  is  steel. 

2)  Cr  radiation  is  used  in  the  analysis.  The  ^-goniometer  geometry 
is  assumed. 

3)  Elastic  constants  do  not  change  with  depth. 

4)  The  constants  a^,  n  _  ,  defined  by  equation  (15)  are  constant 
over  the  maximum  depth  of  penetration  (5.4^). 

5)  For  simplicity,  it  was  assumed  that  the  stress  state  in  the  plane 
of  the  surface  is  isotropic  (c^  *  °22^  and  t^e  variation  Q  n  anc3  ^22 
with  depth  are  identical: 


a°ll_*p21 
dz  bz 


Here,  (cn)  n  was  chosen  arbitrarily  as  -400  MPa. 
X  X  2*U 


6)  Different  exponents  (n_)  were  used  to  form  the  (c,-)  profiles 

33  33 

with  depth.  However,  in  all  cases  values  were  selected  such  that 

y\ 

(c  )  (Eqns-  24,26)  was  approximately  100  MPa.  This  value  was  chosen 
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since  experimental  values  reported  in  literature  on  steels  are  generally 
between  80  -  140  MPa. 

7)  For  a  given  stress  tensor,  the  vs  sin2*  array  was  examined 

as  follows: 


a.  Bi-axial  analysis  (Eqn.  6)  over  the  entire  *-range, 

*«(0P,  6CP),  and,  only  over  the  high  *  Tange,  (45,60). 

b.  Two  tilt*  analysis  using  *  m  0  and  4^,  and  0  and  60°,  and 
45°  and  bCP  respectively. 

c.  Tri-axial  analysis  (Eqn.  24)  over  the  entire  *-range. 

d.  Tri-axial  analysis  over  *  •  (0,  33.21°);  first  four 

points. 

e.  Tri-axial  analysis  over  *  *  (39.23  -  6CP ) ;  last  four 


points. 

8)  The  results  are  shown  in  tables  (1-3)  which  show  the  *-range 
employed,  the  stress  exponents  (n^)  assumed  to  exist  in  the  penetra¬ 
tion  volume,  type  of  analysis  used  to  obtain  the  stresses  the  error 

resulting  from  the  analysis,  the  type  of  error  and  the  correlation  coeffi¬ 
cient  of  the  regression  line  fitted  to  the  ”d”  vs  sin2*  data. 


Group  I 

For  this  group,  a 22  were  assumed  to  be  constant  over  the  depth 

of  penetration.  Three  types  of  gradients  were  assumed  for  o with  the 
stress  exponent  n^  *  1,2,3  respectively.  The  profiles  with  depth 

(Eqn.  19)  for  each  case  are  shown  in  figure  3a.  For  all  three  cases  the 

/S 

average  value  (O33)  (Eqns.  24,26)  is  approximately  100  MPa,  although  the 
gradient  is  much  steeper  for  n^“  3  compared  to  the  gradient  for  n^  *  1- 

The  stress  tensor  sampled  by  the  x-ray  beam  at  any  *-tilt  will  be  : 


(28) 


*In  ’’two-tilt”  analysis,  it  is  assumed  that  all  ”d”  vs  sin2*  values  fall 
on  a  straight  line.-  Thus  only  the  first  and  last  data  points  are  used  to 
obtain  the  stress.  ' 


The  "d"  vs  sinaJr  plots  corresponding  to  stress  tensor  (28)  for  different 
n^  values  are  shown  in  figure  3b,  It  is  seen  that; 

i)  Curvature  in  "d”  vs  sinajr  increases  with  increasing  exponent  n^. 
ii)  All  three  curves  converge  as  sinaJr  increases , (penetration  depth 
decreases)  as  the  gradients  in  stress  also  converge. 

Analysis  for  all  three  cases  (n^  *  1*2,3)  are  summarized  in  table  I. 
Examination  of  this  table  reveals  that: 

1)  If  a  bi-axial  assumption*  is  used  in  analyzing  "d"  vs  sin2jr  data 
obtained  from  surface  layers,  when  a  tri-axial  stress  tensor  is  present, 
appreciable  error  will  occur  in  the  stress  value  calculated  from  the  slope. 

The  error  will  be  the  sum  of  two  components: 

/\ 

1)  The  (S33)  value  that  will  contribute  directly  to  the  slope  (Eqn.  24). 
ii)  The  curvature  error  which  arises  from  fitting  a  least-squares  line 

to  curved  data. 

2)  If  a  tri-axial  solution  is  used  (Eqn.  24),  the  p.-  calculated  from 

A  11 

the  slope  will  be  corrected  for  (0^3) >  but  the  curvature  error  remains. 

3)  Curvature  error  increases  with  increasing  n^,  as  predicted  by  eqn.  25  and 
may  be  quite  high  if  a  low  f  range  is  used. 

4)  The  curvature  error  is  minimized  if  a  least-squares  line  is  fitted 
to  high  Jr  points  only.  This  effect  is  due  to  decreasing  absolute  error, 
associated  with  high  Jr  points,  as  predicted  by  equation  (25). 

It  can  be  seen  from  equation(25)  that,  when  no  gradients  are  present 
in  <7^  and  ,  curvature  error  is  only  a  function  of  the  033  gradient. 

It  is  not  dependent  on  the  magnitude  of  p...  In  figure  4,  the  variation 

/n  11 

of  curvature  error  with  magnitude  of  (O33)  is  plotted  for  n^  *  1,2,3  re¬ 
spectively,  when  a  high  Jr  range  (Jre 39. 23-6CP )  is  used.  Also  included,  as 

a  basis  for  comparison,  is  the  total  error  in  the  stress  obtained  from  a 

A 

bi-axial  analysis  using  the  total  Jr-range  (0,  60P )  for  the  same  (033) 
values  (n^  "  1  for  this  case).  It  is  seen  that,  using  tri-axial  analysis 
over  the  high  Jr -range  will  result  in  a  maximum  curvature  error  of  25  MPa 
(for  (033)  <  170  MBa)  which  is  generally  within  the  total  experimental  error. 
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Groups  II  and  III 

To  check  the  generality  of  the  conclusions  arrived  at  in  group  I,  the 
effects  of  stress  gradients  in  other  terms  of  the  stress  tensor  were  inves¬ 
tigated: 


In  group  II,  a  stress  tensor  of  the  form, 


(c 


ijN 


(29) 


was  assumed,  where  (a m  -400  +  K^z2  MPa, 


The  calculated  stress  profile  with  depth  for  this  case  is  shown  in 
figure  5a,  with  the  resulting  Mdn  vs  sin2ty  in  figure  5b.  The  curvature 
agrees  with  that  calculated  by  Dolle  et  al.  However  high  $  values  are 
needed  to  see  the  effect,  since  the  gradient  is  not  very  steep  compared 
to  those  in  figure  3a. 

w. 

Results  of  analysis  are  shown  in  table  II.  The  (9^3)  values  obtained 
in  this  case  are  residuals  due  to  misfit  of  the  linear  regression  equation. 
When  a  high  i|r- range  is  used,  the  residual  is  well  within  experimental  error. 
It  is  seen  that,  whereas  analysis  of  tri-axial  data  by  bi-axial  methods 
yields  very  high  error,  analysis  of  bi-axial  data  by  tri-axial  methods  is 
more  accurate,  especially  if  the  high  i(r-range  is  employed. 

In  group  III  the  effect  of  combined  gradients  were  examined.  Gradients 
were  assumed  in  all  components  of  the  stress  tensor  (13).  Two  cases  were 
examined: 


a)  <Cn>^  "  -400  +  Kj^z,  where  K,,  -  5  MPa/p, 


'll 


"  K33z8 


b>  <°11>*  "  -400  +  Knz2 


<°33>*  "  K33z3 
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The  calculated  stress  profiles  for  both  cases  are  shown  in  figures  (6a) 
and  (6b)  respectively.  The  resultant  "d"  vs  sin2i|r  graphs  are  shown  in  figure 
(6c),  Again  both  lines  converge  for  high  f  angles.  The  results  of  analysis 
are  summarized  in  table  III, 

Examination  of  tables  II  and  III  show  that  the  conclusions  from  group 
I  are  indeed  general.  It  can  also  be  seen  that  a  high  correlation  coeffi¬ 
cient  in  MdM  vs  sin3f  is  not  a  good  measure  of  the  dimensionality  of  the 
stress  tensor. 

Summary 


1)  Deviations  of  interplanar  spacing  "d"  vs  sin2^  from  linearity, 
arising  from  stress  gradients  are  generally  small,  confirming  previous 
comments  on  this  matter. ^ 

2)  Analysis  of  even  slightly  curved  "dn  vs  sin2ijf  plots  assuming  a 
bi-axial  surface  stress  state  may  cause  high  error  in  the  stress  value 
determined  from  the  slope, 

3)  Employing  the  methods  for  tri-axial  stress  tensors  described  here 
will  minimize  errors  due  to  the  average  value  of  the  stress  normal  to  the 
surface.  There  will,  however,  still  be  some  residual  error  associated  with 
curvature  of  the  data, 

4)  This  curvature  error  may  be  minimized  by  using  high  f  points  in 
the  linear  least-squares  fit,  (sin2^  >,4). 

5)  When  two  radiatfons  are  available  for  a  measurement,  the  one  having 

a  high  absorbtion  coefficient  (low  penetration  depth)  should  be  used  if  stress 
gradients  with  depth  are  present.  This  will  also  minimize  errors  due  to 
curvature. 

6)  Unless  it  is  absolutely  clear  that  stresses  normal  to  the  surface 
are  absent,  or  quite  small,  the  current  x-ray  method  of  stress  analysis 
should  be  replaced  by  the  methods  described  here. 
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TABLE  TITLES 


Table  I:  Analysis  of  ,,dM  vs  sin3*  plots  1-1,  1-2,  1-3,  in  figure  3b.* 

Table  II:  Analysis  of  "d"  vs  sin3*  shown  in  figure  5b.** 

Table  III:  Analysis  of  "d"  vs  sin3*  plots  III-l,  III-2,  shown  in  figure  6b.*** 


*  For  ell  cases  is  -400  MPa  <a^>  values  shown  are  average  values; 

is  zero  at  Hie  surface  z-0,  Tne  gradient  in  (o ^ 3>,  is  shown  in 
figure  3a  for  all  three  cases,  n^*l,2,3. 

**  -  a^AzM})  is  -400  MPa  and  g33(z)-0.  The  (o33>  values  obtained 

from  the  analysis  are  residual  errors  due  to  the  misfit  of  the  linear 
model  used  in  analysis. 

***In  this  case  there  are  gradients  in  both  c1T  and  a..,  (figure  6a),  a  ° 
-  O11(z-0)  -  -400  MPa  and  a33(z«0)  -  0.  11 


TABLE  I 


TABLE  II 


* -Range 
(Degrees) 

nll 

Type  of 
Analysis 

<~n> 

(MPa) 

.  A 

<033> 

(MPa) 

Error  in 

i<5i>  ■  ffu°i 

Type  of 

Erron  in 

Correlation 
Coefficient  of 
"d"  vs  sin2* 
plot. 

0-45 

n 

Si-axial 

"Two-Tilt" 

-336.4 

0 

63.6 

Curvature 

0-60 

n 

II 

-354.5 

0 

45.5 

»t 

45.60 

n 

II 

-361.8 

0 

38.2 

n 

0-60 

B 

Bi-axial 

"sin2*1' 

-355.66 

0 

44.34 

ii 

-.9993 

45-60 

n 

„ 

_ i 

-415.8 

0 

15.8 

it 

-.9987 

0-60 

B 

252311 

-364.86 

-9.20 

35.14 

it 

-.9993 

0-33.21 

B 

ti 

-305.60 

24.53 

94.40 

ii 

-.9999 

39.23-60 

B 

n 

-382.61 

-2.272 

17.39 

ii 

-.9986 

TABLE  III 


t -range 
(Degrees) 

nll 

°33 

Type  of 
Analysis 

<^i> 

(MPa) 

<C33> 

(MPa) 

Eryor  in 

|f<air>1“®Li°i 

Type  of 
error  in 

•a  ! 

Correlation 
Coefficient  of 
•d"  vs  s  in2  $  plot. 

0-45 

D 

Bl-axial 

MTwo-TiltM 

-575.5 

0 

175.5 

m 

0-60 

B 

2 

II 

-535.7 

0 

135.7 

ii 

45-60 

B 

2 

II 

-456.0 

0 

56.0 

ii 

0-60 

B 

2 

Si-axial 

"sin3*" 

-540.08 

0 

140.08 

-.9982 

45-60 

B 

2 

II 

-458.4 

0 

58.4 

ti 

-.9996 

0-60 

B 

2 

Tri-axial 

-436.90 

103.20 

36.90 

Curvature 

-.99o2 

0-33.21 

B 

2 

-583.70 

19.57 

183.70 

it 

-.9999 

39. 23-60 

B 

2 

ii 

-393.39 

86.36 

6.61 

i» 

-.9992 

0-45 

2 

2 

Bi-axlal 

"Two-Tilt" 

-586.1 

0 

186.1 

MSKM 

1353533533 

0-60 

2 

2 

II 

-550.0 

0 

150.0 

n 

45-60 

2 

2 

II 

-477.6 

0 

77.6 

n 

0-60 

2 

2 

Bl-axial 

"sin3*" 

-548.29 

0 

148.29 

n 

-.9985 

45-60 

2 

2 

II 

-477.3 

0 

77.3 

n 

-.9999 

0-60 

2 

2 

-431.98 

116.31 

31.98 

Curvature 

-.9985 

0-33.21 

2 

2 

n 

-565.30 

40.42 

165.30 

ii 

-.9997 

39.23-60 

2 

2 

n 

-392.02 

100.73 

7.98 

n 

-.9995 
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FIGURE  CAPTIONS 


Figure  1  : 

Figure  2  : 

Figure  3a: 

Figure  3b: 

Figure  4  ; 


Figure  5a: 

Figure  5b: 

Figure  6a: 
Figure  6b: 
Figure  6c: 


Definition  of  the  angles  $  and  |  and  orientation  of  the  labora¬ 
tory  system  with  respect  to  the  sample  system  and  the 
measurement  direction  S.» 

Penetration  depth  ,,Tn  vs  sin2^  for  Cr-radiation  on  steel  at 
20  -  156°. 

Calculated  profiles  vs  depth  for  group  I.  In  all  three 

cases  0^  m  0 22  "  MPa  throughout  the  depth  of  penetration. 

The  !,dn  vs  sin3ijr  corresponding  to  the  stress  profiles  shown  in  figure 
3a.  (The  results  of  analysis  for  these  curves  are  shown  in  Table 
I,  groups  1-1,  1-2,  1-3  respectively.) 

Variation  of  the  magnitude  of  error  in  the  surface  stress 
with  the  magnitude  of  the  average  stress  in  the  direction  of 
the  surface  normal.  The  steepest  curve  shows  the  error  in  sur¬ 
face  stress  determined  by  bi-axial  analysis  for  <|r  *  0-6CP  .  The 
other  three. curves  show  errors  resulting  from  triaxial  analysis 
for  t  -  39.23  -  60 P. 

Residual-stress  profile  vs  depth  for  group  II.  In  this  case 
there  is  no  stress  gradient  in  the  direction  of  the  surface 
normal. 

3 

Xhe  "d"  vs  sin  f  plot  corresponding  to  the  stress  profile  with 
depth  shown  in  figure  5a.  Results  of  analysis  of  this  curve  are 
summarized  in  Table  II. 

Residual  stress  profiles  with  depth  for  group  III-l. 

Residual  stress  profiles  with  depth  for  group  III-2. 

The  "d11  vs  sin3^f  corresponding  to  residual  stress  profiles  shown 
in  figures  6a  and  6b.  The  results  of  analysis  for  both  curves 
are  summarized  in  Table  III,  groups  III-l,  III-2  respectively. 
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